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Introduction 


FINAL  REPORT 


The  association  of  intensity  variation  in  atmospheric 
emissions  with  the  presence  of  atmospheric  waves  has  been  made 
as  far  back  as  1957  by  Krassovsky  (Krassovsky  (1957)  ).  Later, 
these  effects  have  been  observed  and  analyzed  by  a  number  of  in¬ 
vestigators  Qe.g.  Okuda  (  1962  ),  Silverman  (  1962  ),  Barbier  (196’: 
1965),  Weill  and  Christophe-Glaume  (1967),  Dachs,  (1968),  Andrews 
(1976),  Dyson  and  Hopgood  (1978  )j|»  for  both  the  6300  8  01  and  the 
5200  8  NI ,  The  presence  of  gravity  waves  in  OH  emissions  has  been 
examined  by  Krassovksy  (1972),  Krassovsky  and  Shagaev  (1974), 
Krassovsky  et  al  (1975),  Armstrong  (1975),  Krassovsky  and  Shagaev 
(1977),  Moreels  and  Herse  (1977),  and  Peterson  (1979).  Observa¬ 
tions  'm  the  effect  of  gravity  waves  on  02(^5:)  emissions  have  been 
made  by  Noxon  (1978) .  Theoretical  treatments  for  the  effects  of 
gravity  waves  on  airglow  have  been  made  by  Porter  et  al  (1974)  for 
6300  8  01  and  5200  A  NI  and  Weinstock  (1978)  for  0 2 ( 1 Z )  and  OH 
emissions . 

In  general,  the  treatment  of  gravity  wave  effects  on  air- 
glow  require  an  analysis  of  the  effects  of  gravity  waves  on  the  in¬ 
dividual  atmospheric  constituents  which  produce  the  airglow  through 
chemical  reactions.  There  are  two  cases  to  be  considered.  The 
first  is  when  the  atmospheric  constituents  arc  minor  components  of 
the  atmosphere.  It  is  then  usual  to  assume  that  the  gravity  waves 
are  carried  by  the  major  atmospheric  constituents  and  their  effects 
on  each  of  the  minor  const i  tuo-'  ts  are  analyzed  with  the  assumption 
that  the  effects  of  the  minor  constituent  on  the  gravity  waves 
are  negligible.  The  second  case  is  when  the  atmospheric 
constituent  producing  the  particular  airglow  emission  is  actually 
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itself  the  major  constituent.  In  this  case  there  is  no  longer  a 
need  to  separately  calculate  the  response  of  the  constituent  to 
the  gravity  waves.  The  problem  of  calculating  airglow  is  then 
made  much  simpler  provided  a  suitable  gravity-wave  model  is  chosen. 

In  calculating  the  effects  of  gravity  waves  on  the  oxygen 
red  line  and  the  nitrogen  green  line,  it  is  necessary  to  first 
consider  the  effects  of  gravity  waves  on  charged  particles  such  as 
the  electron,  since  it  is  primarily  through  dissociative  recombina¬ 
tion  of  these  electrons  with  positive  ions  (C^*,  N0+)  which  pro¬ 
duce  these  particular  airglow  emissions.  The  problem  is  further 
complicated  by  the  presence  of  the  magnetic  field  which  puts  ad¬ 
ditional  constraints  on  the  motion  of  the  charged  particles. 

Detailed  knowledge  of  the  response  of  the  charged  particles  to  the 
gravity  wave  including  phase  relations,  magnitude  of  response  etc. 
have  been  treated  by  Thome  (1968),  Testud  and  Francois  (1971), 
Klostermeyer  (  1972a, b)  and  Porter  and  Tuan  (1974).  With  the  ex¬ 
ception  of  Thome  (1968),  all  the  rest  make  use  of  gravity-wave 
models  which  include  dissipation  and  also  take  into  consideration 
the  diffusion  (ambipolar)  of  the  charged  particles.  The  calculations 
include:  (1)  the  magnitude  of  the  response  which  depends  on  the 

vertical  gradient  of  the  undisturbed  ionospheric  number  density 
profile  for  the  minor  constituent,  (2)  the  phase  variation  with 
height  of  the  response  relative  to  the  gravity  waves. 

Mere  we  should  mention  that  it  was  Thome  (1968)  who  first 
observed  a  180°  phase  change  in  the  response  of  the  F-layer  above 

t 

and  below  the  F-region  peak  and  provided  a  suitable  physical 
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explanation.  Subsequent  theoretical  papers  by  Testud  and  Francois 
(1971),  Klostermeyer  (1971a, b)  and  Porter  and  Tuan  (  1974  )  have 
quantitatively  confirmed  the  phase  variation. 

To  calculate  the  effects  of  gravity  waves  on  OH  emission 
it  is  necessary  to  know  the  effect  of  gravity  waves  on  neutral 
minor  atmospheric  constituents  such  as  H  and  0^.  The  effects  of 
gravity  waves  on  neutral  atmospheric  constituents  have  been  con¬ 
sidered  by  Dudis  and  Reber  (1976)  who  use  the  Hines  isothermal 
gravity-wave  model  and  neglect  wave-induced  diffusion  to  obtain 
simple  analytic  phase  relationships  between  the  response  of  the 
atmospheric  constituent  and  the  gravity  wave.  The  method  is  just¬ 
ifiable  for  lower  atmospheres  where  dissapation  is  not  crucial  and 
for  free  gravity  wave  modes  with  low  horizontal  phase  velocity 
(less  than  100  m/sec.).  Chiu  and  Ching  (1978)  have  considered  the 
effect  of  gravity  waves  on  minor  constituents  with  a  layered 
Structure.  Like  Dudis  and  Reber  (1976)  they  have  used  the  Hines 
(1960)  analytic  model  for  the  gravity  wave  and  obtained  similar 
relationships  in  terms  of  well-known  atmospheric  parameters  such 
as  the  scale  height  and  the  ratio  of  specific  heats,  etc.  While 
their  conclusions  are  essentially  similar  to  those  mentioned  in 
the  previous  paragraph,  they  also  specifically  point:  out  that  the 
response  for  the  bottom  side  is  always  greater  than  the  top  side. 
This  last  result  is  only  strictly  true  for  a  Hines'  gravity  wave 
model  in  which  the  rate  of  increase  of  the  horizontal  and  vertical 
velocity  components  with  height  are  always  the  same. 

In  this  paper  we  wish  to  consider  the  effect  of  gravity 
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waves  on  OH  emissions  using  a  gravity-wave  model  JjTuan  (1976)] 
which  satisfies  the  rigid  surface  boundary  condition  at  the 
ground  level  and  can  provide  for  guided  as  well  as  free  modes. 

The  undisturbed  hydrogen  and  ozone  profiles  are  taken  from  Good 
(1976)  and  Keneshea  and  Zimmerman  (private  communication).  We 
will  consider  the  effect  of  gravity  waves  on  the  OH  emission 
profile  for  different  horizontal  phase  velocities.  It  will  be 
shown  that  for  sufficiently  low  horizontal  phase  velocity  (~  33 
m/scc.)  the  gravity  wave  can  produce  structures  in  the  OH  emission 
profile  and  that  the  structures  tend  to  but  need  not  always  be  more 
pronounced  on  the  bottom  side.  These  structures  reveal  the 
presence  of  dark  areas  usually  observed  below  the  peak  intensity 
which  appear  to  correspond  to  the  "holes"  observed  by  Peterson 
(1978).  We  will  consider  the  "phase"  relationship  between  the  re¬ 
sponse  of  H  and  0^  and  the  gravity  wave.  In  this  case,  just  as 
for  the  case  of  the  ionosphere,  the  "phase"  varies  continually 
from  the  bottom  side  to  the  top  side.  Here,  we  use  the  word  "phase" 
only  in  a  very  approximate  senbe,  since  neither  our  gravity-wave 
model  nor  the  response  of  the  atmospheric  constituents  have  a 
simple  periodic  spatial  variation. 

For  the  case  of  large  horizontal  phase  velocities 
(say,  >  160  m/scc.),  the  vertical  "wave  length"  is  then  too  large 
compared  wi th  the  thickness  of  the  ozone  layer  to  produce  any 
significant  spatial  structures.  However,  there  will,  of  course,  be 
significant  temporal  variation  in  the  total  columnar  intensity 
wti  i  ch  can  be  observed  by  a  ground-based  photometer. 
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Theoretical  Formulation 


We  shall  assume  that  the  primary  mechanism  for  the  pro¬ 
duction  of  OH  emissions  is  given  by  the  reaction 

H  +  03  OH  (v)  +  02  (1) 

The  OH  intensity  profile  I  (photon/km^sec)  is  then  given  by: 

I  "  )  21k(v)  [5I]  [03J  (2) 

where  J  is  the  efficiency  of  photon  emission  per  Oil 
molecule  formed  and  k(v)  is  the  rate  coefficient  for 
equation  (1).  It  is  possible  to  show  that  in  regions 
between  80  km  and  100  km  [Earl  Good  (1976)  ,  Zimmerman 
(private  communication)]  chemical  equilibrium  is  a  good 
assumption  for  the  production  and  loss  of  ozone.  Thus, 
from  equation  (2),  all  we  need  is  to  determine  the  effect 
of  gravity  waves  on  the  hydrogen  and  ozone  profiles.  Since 
both  are  obviously  very  much  minor  atmospheric  constituents, 
corresponding  to  the  first  case  mentioned  in  the  Introduction, 
we  shall  assume  that  the  gravity  wave  is  being  carried  by  the 
major  constituents  and  consider  its  effect  on  H  and  O3 
separately.  We  shall  also,  for  the  present  paper,  make  use 
of  the  linearized  (perturbation)  approach  and  will  specify 
other  approximations  as  they  appear.  The  continuity  equa¬ 
tion,  momentum  conservation  equation  for  any  given  minor 


atv.io.-.ph.- r •  :n.  lows : 

-  'I  ~  ^  —  eL-\r  CN^r) 

U  t  (3) 

mN>?  4  ^N^CV-AT?)  •=  -fvf-xvx/vrl 

J)t  t  l  0  JC4) 

concentration  of  the  minor  component 
velocity  field 

molecular  mass  of  the  component 

collision  frequency  between  the  component 
and  the  major  atmospheric  constituents 

partial  pressure  of  the  minor  component 

gravity-wave  velocity  field  given  by  the 
gravity-wave  model 

production  rate 

loss  rate 

Even  with  the  assumption  of  a  model  for  ,  equations  (3) 

0 

and  (4)  are  not  closed.  An  additional  condition  for  the 
equation  of  state  is  needed  to  close  the  system  of  equations. 
However,  in  general  we  may  need  different  equations  of  state 
for  the  perturbed  and  the  unperturbed  atmospheric  constituent 
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For  instance,  if  the  perturbation  is  an  oscillation,  we  may 
require  the  adiabatic  condition  which,  however,  cannot  be 
validly  applied  to  the  unperturbed  atmosphere.  Thus,  we 
will  leave  the  equations  open  and  close  them  with  one  or 
more  appropriate  equations  of  state  as  the  need  arises. 

We  will  assume 


N  =  N0  + 

t  ■  V 

L  -  L.  ♦  at 

t  ■  h  *  Ar 

r  -  n  *  ^ 


(5) 


where  NQ ,  ,  and  are  the  background  concentra¬ 

tion,  etc.  of  the  minor  constituent  undisturbed  by  the  gravity 
wave.  ~VT^  is  the  drift  velocity  of  the  minor  constituent 
through  the  background  atmosphere  and  is  always  strong  whenever 
there  is  a  sharp  concentration  gradient.  We  shall  assume  that 
these  dependent  variables  are  independent  of  time  and  are 
functions  of  only  (horizontal  stratification).  A  N ,  A Jj.  , 
AL  ,  ,  and  are  the  perturbed  dependent 

variables  produced  by  the  gravity  wave.  We  shall  also  neglect 
variation  in  collision  frequency. 

Substituting  equation  (5)  into  (3)  and  (4)  we  obtain  the 
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following  zeroth  order  equations 
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horizontal  wave  vector.  The  equations  become, 
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We  shall  now  consider  the  right-hand  side  of  equations  (9) 
and  (10).  The  collision  frequency  p  at  100  km  is  of  the 
order  of  2  x  10^  sec'^,  while  the  gravity  wave  angular 
frequency  for  a  2  hour  period  is  of  the  order  of  8.7  x  10"^ 
sec'*.  Hence, 

4  x  10~7  (11) 

P  o 


We  can  immediately  neglect  the  two  velocity  terms, 
term  involving  q  ,  we  have 
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and  AW<  is  of  the  order  of  0.5  m.  sec.’ 
so  we  can  neglect  this  term. 

The  pressure  term  in  equation  (9)  may  be  written  as 


Now 


AN 

7v, 


<  1 


(13) 


where  j)  is  the  diffusion  coefficient  of  the  gas. 
The  drift  velocity  term  is  given  by 


X>  (AH  \ 
Nr  A  h  [  NcJ 


(  H) 


In  general,  the  diffusion  coefficient  T)  is  of  the  order  of 

200  m^sec  *  at  100  km.  For  the  bottom  side  of  ozone  •  r~^~  - 

Nt  J  f- 

0.4  x  10  ^  m ,  *  Hence  3  x  10  ^  m.sec  * . 


Since 


^  i  i 

i-f  '  1  ’ 


we  may  neglect  this  term  in  comparison  with 
The  pressure  term  in  equation  (9)  may  be  written  as 


h 


(15) 


Since  maximum  k  <v->  6  x  10 
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and 
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m.sec 


which  is  again  negligible. 
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Li 


The  term 
written  as 


_J _  W) 

**-N6  Jc  l  t- 


^  N<  X  6  h 


can  be  approximately 


where  we  assume  that  to  a  rough  approximation,  can 

be  approximated  by  a  sinusoidal  wave  with  vertical  wave- 
vector  k,  •  The  maximum  k±  we  consider  correspondinj 


vector  k ,  .  The  maximum  k±  we  consider  corresponding 

to  a  horizontal  phase  velocity  of  about  33  m/sec,  is  of  the 
order  of  k±  ~  1.3  x  10'3  m.“*.  This  means  that 

0.26  m.sec.  * ,  which  is  again  very 

much  less  than  the  typical  value  of  4W1  of  some  m.sec.  * . 
Thus,  for  all  horizontal  phase  velocity  of  interest,  we  may 
assume  that  equations  (9)  and  (10)  reduce  to 

Ate.  —  Au.  I 

AvJ  - 


Substituting  equations  (17)  and  (18)  into  equation  (8)  and 
assuming  that  A L  we  obtain: 

an  ^  ^  -  i-  -l 

CO  <J  °°  <)  & 


M  Au.  +-  — 
)  ^  ^ 


Thus,  given  any  gravity-wave  model,  we  can  calculate  Ah] 
from  equation  (19) . 
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Gravity-wave  model 

The  gravity-wave  model  we  use  [Tuan  (1976)]  satisfies 
the  rigid  surface  boundary  condition  at  the  ground  level 
and  propagates  in  an  atmosphere  in  which  the  variations  in 
temperature,  mean  molecular  mass,  acceleration  due  to  gravity 
and  ratio  of  specific  heats  are  all  taken  into  account  using 
the  1972  COSPAR  model  wherever  feasible.  The  model  also 
allows  for  the  propagation  of  discrete  as  well  as  for  the 
continuous  or  free  modes.  Since  we  are  only  concerned  with 
altitudes  below  110  km,  dissipation  due  to  thermal  conduction, 
viscosity  and  ion  drag  may  be  neglected.  In  the  complex 
notation,  ^  is  actually  pure  imaginary,  so  &N 
is  actually  real  inequation  (19).  In  general,  the 
boundary  condition  at  the  ground  ensures  that  the  dependent 
variables  (i.e.,  or  )  »  which  are  functions  of  ir 

only,  are  either  pure  real  or  pure  imaginary.  Since  we  use 
a  complex  time  dependence  [exp  (  i  wi  )]  ,  this  means  that 
there  is  a  stationary  wave  along  the  vertical  direction 
in  contrast  to  the  travelling  wave  along  the  horizontal 
direction,  an  expected  result  for  a  rigid  surface  boundary 
condition. 

The  model  we  use  [Tuan  (1976)]  actually  only  calculates 
the  pressure  variation  given  by 

ifs  -  ^k 

I  —  - 1~  (20] 

w* 

where  J*o  is  the  background  density  and  4  j*  is  the  pressure 


variation  of  the  gravity  wave.  The  velocity  field  may  be 
obtained  from  the  pressure  field  through  the  hydrodynamic 
equations  and  the  results  are  given  by 

p  -  4  * 

aw>  -  - Morl 

where  Au>  -  (J2  t-sL 

1  ^c.*-  c  (j-  J 

In  equation  (21)  ,  since  ^  obeys  the  boundary  condition 
[Tuan  (1976)] 


-  M°)  f 


The  vertical  velocity  field  obeys  the  usual  well-known 
boundary  condition  of  AWCe)  =  0. 

Figs  (J.)  and  (2)  show  the  vertical  and  horizontal 
velocity  fields  A  U/  and  .  A  U  as  functions  of  A  ,  the 
altitude,  for  horizontal  phase  velocities  of  302 
and  33  m/sec,  respectively.  For  both  the  vertical  and 
horizontal  velocity  fields  there  are  nodes  along  the 
vertical  direction.  This  moans  that  there  are  adjacent 
layers  where  the  particles  move  in  opposing  directions. 

The  approximate  vertical  wave  lengths  are  of  the  order  of 
150  -  200  km  and  8  km,  respectively.  Since  the  tot3l  half 


width  of,  for  instance,  the  hydrogen  number  density  profile 
in  of  the  order  of  10  km,  spatial  structures  arc  expected  only 
for  the  gravity  w.ive  with  a  horizontal  phase  velocity  of  33  m/sec 
(or  lens).  Here,  we  should  mention  that  the  302  m/sec .  ^  wave 
is  very  close  to  the  guided  Lamb  mode  so  that  "vertical  wave 
length"  is  meaningless  below  100  km  altitude,  where  the  length 
scale  in  the  atmospheric  structure  is  less  than  the  wave  length 
and  the  Hines  (I960)  model  is  no  longer  valid.  At  least  for  this 
altitude  range  (80  -  105  km)  the  horizontal  velocity  field 
increases,  with  height  far  more  rapidly  than  the  vertical  velocity 
field.  In  Tig.  (1)  our  gravity-wave  model  shows  a  maximum 
hor'i  zontnl  velocity  of  84  m/sec."1  at  91  km.  This  is  consistent 
with  Peterson’s  estimation  (1979)  of  71  m/sec."^  at  an  assumed 
90  km  altitude  from  an  observed  OH  emission. 


Result  s 

For  the  unperturbed  hydrogen  and  ozone  number  density 
profile's,  we  use  the  semi-empirical  data  of  Good  (  1976  )  and 
KenesluM  and  Zimmerman  (private  communication).  In  Figs.  (3)  and 
(4)  we  show  l he  "unperturbed"  variations  in  the  ozone  and  hydrogen 
number  density  as  a  function  of  altitude.  The  Zimmerman  data 
is.  used  throughout  for  both  Figures.  We  sec:  that  the  hydrogen 
concentration  has  a  peak  at  between  84  to  85  km  while  the  ozone 
has  a  local  peal  at  about  86  km. 
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The  concentration  gradient,  at  least  for  the  hydrogen 
(Fig.  5)  curve,  is  steeper  below  the  peak  than  above  the 
peak . 

Fig.  (5)  shows  variations  in  ozone  AhLO^^ 
hydrogen  AN  (J"0  and  total  ANj.  number  densities  when 
a  gravity  wave  with  a  horizontal  phase  velocity  of  302  m.sec'1 
and  a  two-hour  period  passes  by.  In  Figs.  (5)  to  (10) 
(inclusive)  the  number  density  variations  AN  are  all 
computed  from  equation  (19).  For  Figs.  (5)  to  (9)  (inclusive) 
we  have  used  Zimmerman's  data  for  the  unperturbed  concen¬ 
tration  profile  N0.  For  Fig. (10)  we  have  used  Earl  Good's 
(1976)  data  for  the  unperturbed  profile.  The  two  terms  on 
the  right-hand  side  of  equation  (19)  represent  the  individual 
contributions  from  the  horizontal  and  vertical  diffusion 
to  AN  .  For  both  the  minor  constituents  [5(a)  and  5(b)] 
the  vertical  diffusion  is  predominant  below  the  peaks.  Since 
the  horizontal  diffusion  has  the  same  sign  (negative)  as 
the  vertical,  there  is  significant  reinforcement.  The 
vertical  diffusion  tends  to  dominate  below  the  peak  because 
of  the  relatively  steep  vertical  concentration  gradient.  At 
the  same  time,  the  horizontal  velocity  field  below  the  peaks 
is  not  yet  sufficiently  greater  than  the  vertical  that 
horizontal  diffusion  can  effectively  compete  with  the 
vertical  diffusion.  We  might  mention  that  because  the 
horizontal  and  vertical  diffusion  reinforce  each  other  for 
both  minor  constituents,  the  response  is  significantly 
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greater  below  peak  for  "V7l„  =  302  m.sec"*  than  for  other 
gravity  waves  tvith  lower  horizontal  phase  velocities. 

For  £h/r  there  is  no  sharp  vertical  concentration  gra¬ 
dient  (no  peaks)  for  the  unperturbed  total  atmospheric 
concentration  and  the  horizontal  diffusion  dominates. 

Unlike  the  minor  constituents,  there  is  considerable 
cancellation  between  vertical  and  horizontal  diffusion  in 
A/^-p  ,  which  also  accounts  in  part  for  the  very  large 
response  (i.e.,  large  ratios  for  AAJCOj) 

and  A  N  CH) y/  ANp-  as  compared  with  other  gravity  waves. 

Above  the  hydrogen  and  ozone  peaks  vertical  domination 
continues  for  ozone  but  for  hydrogen,  which  has  a  less 
steep  vertical  gradient,  horizontal  diffusion  eventually 
becomes  dominant.  The  ANj-  of  course  continues  to  have 
large  cancellation  with  a  horizontal  domination. 

Fig.  (6)  shows  AN(03),  AN(H)  and  A  NT  for  a  gravity  wav 
with  a  horizontal  phase  velocity  of  78  m.sec.'*  Below 
the  peaks,  we  again  have  high  vertical  domination  for  the 
minor  constituents  with  a  horizontal  domination  at  lower 
altitudes  for  AN^.  Here  cancellation  occurs  to  a  much 
less  extent  than  for  the  previous  case.  Above  the  peaks, 
horizontal  domination  eventually  takes  over  for  both 
ozone  and  hydrogen.  In  both  Fig.  (6a)  and  (6b),  one  can 
begin  to  see  a  phase  change  in  the  response  of  the  minor 
constituent  relative  to  An^..  Below  the  peaks  the 
response  is  out  of  phase,  but  above  the  peaks  they  appear 
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and  ANOO/^Nr  as  compared  with  other  gravity  waves. 

Above  the  hydrogen  and  ozone  peaks  vertical  domination 
continues  for  ozone  but  for  hydrogen,  which  has  a  less 
steep  vertical  gradient,  horizontal  diffusion  eventually 
becomes  dominant.  The  Ah/y  of  course  continues  to  have 
large  cancellation  with  a  horizontal  domination. 

Fig.  (6)  shows  AN(03),  &  N  (H)  and  A  NT  for  a  gravity  wave 
with  a  horizontal  phase  velocity  of  78  m,sec.‘*  Below 
the  peaks,  we  again  have  high  vertical  domination  for  the 
minor  constituents  with  a  horizontal  domination  at  lower 
altitudes  for  AN^,.  Here  cancellation  occurs  to  a  much 
less  extent  than  for  the  previous  case.  Above  the  peaks, 
horizontal  domination  eventually  takes  over  for  both 
ozone  and  hydrogen.  In  both  Fig.  (6a)  and  (6b),  one  can 
begin  to  see  a  phase  change  in  the  response  of  the  minor 
constituent  relative  to  An^,.  Below  the  peaks  the 
response  is  out  of  phase,  but  above  the  peaks  they  appear 
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to  get  in  phase. 

In  Fig.  (7),  the  vertical  diffusion  again  dominates 

below  the  peaks  for  the  minor  constituents.  For  A 
there  is  horizontal  domination  only  at  lower  altitudes. 

Above  the  peaks,  horizontal  domination  and  some  cancellation 
occur  for  both  hydrogen  and  ozone.  Once  again,  one  can  see 
that  the  response  of  the  minor  constituents  is  out  of  phase 
with  the  major  constituents  below  the  peaks  but  gets  in 
phase  above  the  peaks. 

In  Fig.  (8)  and  (9)  we  show  separately  the  response  of 
and  h  N  as  compared  with  ANt.  With  a 

vertical  wave  length  approximately  in  the  range  of  8  -  9  km 
one  can  see  rather  easily  the  response  which  is  about  180° 
out  of  phase  below  say  87  km  but  gets  almost  completely  in 
phase  above  90  km.  The  gravity  wave  used  here  has  a 
horizontal  phase  velocity  of  33  m.sec."*  and  we  have  used 
both  Zimmerman's  data  (Fig.  8  and  9)  and  Carl  Good's  data 
(Fig.  10)  to  show  up  any  significant  difference  in  the 
response  of  the  minor  constituents.  All  in  all  there  is 
relatively  little  qualitative  difference  in  the  basic 
features.  For  both  sets  of  data,  there  is  either  vertical 
domination  or  large  cancellation  (for  hydrogen  and  Carl 
Good's  data)  below  the  peaks.  Above  the  peaks,  there  is 
either  horizontal  domination  and  large  cancellation  (for 
Earl  Good's  data)  or  alternating  vertical  and  horizontal 
domination  (for  Zimmerman's  data). 
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Finally,  we  use  equation  (2)  to  calculate  the  total 
hydrogen  intensity  profile  for  a  gravity  wave  with 

=  33  in. sec.  For  this  case,  as  already  mentioned, 
the  vertical  wave  length  is  less  than  the  half  width  of  the 
unperturbed  OH  profile  and  one  would  expect  structures 
to  occur.  Fig.  (11)  shows  the  OH  intensity  profile,  using 
Zimmerman's  unperturbed  data.  Fig.  (12)  shows  the  profile 
using  Earl  Good's  (1976)  data.  We  see  immediately  for  both 
curves  a  double  layered  structure  below  the  peak  and  a 
suggestion  of  some  kind  of  structure  above  the  peak.  It  is 
obvious  that  the  structure  below  the  peak  is  far  more 
pronounced  than  above  the  peak  in  agreement  with  the  analysis 
of  Chiu  and  Ching  (1978)  in  which  the  Hines  (i960)  gravity- 
wave  model  was  used.  The  "valley"  immediately  below  the 
peak  can  easily  correspond  to  the  "dark  areas"  observed  by 
Peterson  (1978).  While  we  have  found  that  it  is  sometimes 
possible  (depending  on  the  phase  of  the  gravity  wave)  to 
observe  structure  above  the  peak  OH  intensity  rather  than 
below,  it  is  more  likely  to  find  highly  pronounced  structures 
below  peak  intensity. 

In  our  computing  procedure,  the  only  real  difficulty 
comes  from  numerical  differentiation.  In  obtaining  A  N 
in  equation  (19),  we  essentially  have  to  differentiate  the 
function  ^  given  by  equation  (20)  twice,  once  through 
equation  (21)  and  once  through  equation  (19)  .  Various 
smoothing  procedures  have  been  adopted  but  were  found  to 


be  unreliable  because  they  tend  to  alter  the  phase  of  the 
final  ZjM.  In  tn_  nu,  we  have  leciied  to  plot  the 
curves  in  Fig.  5-1M  without  anv  smoothing.  Much  of  the 
structures  in  AN-p  in  Fig.  5(a)  and  Fig.  5(b)  is  caused 
by  small  irregularities  in  the  original  function  'f’  wh  i  c  h 
become  greatly  exaggerated  when  it  is  twice  differentiated. 
These  irregularities  become  much  less  influential  when 
itself  has  a  much  smaller  vertical  wave  length,  as  can  be 
seen  in  Fig.  (8),  (9),  and  (10)  in  which  are 

essentially  smooth  oscillating  curves. 

Conclusion 

The  principal  features  of  the  present  investigation 
may  be  summarized  as  follows:  (1)  We  have  found  that  the 
variations  N  in  number  densities  of  hydrogen  and  ozone 
produced  by  gravity  waves  can  be  attributed  primarily  to 
wave  induced  vertical  diffusion  below  the  layer  peak  and  a 
mixture  of  horizontal  and  vertical  diffusion  above  the  peak. 
For  the  most  part,  the  contribution  of  this  vertical 
diffusion  is  negative  creating  a  decrease  in  the  number 
density  from  the  unperturbed  profile.  The  vertical  dif¬ 
fusion  dominance  can  be  in  part  attributed  to  the 
relatively  steep  vertical  gradient  for  the  minor  constituents. 
The  increase  in  importance  of  the  horizontal  diffusion  above 
the  peak  can  be  in  part  attributed  to  the  good  deal  more  rapid 
increase  with  height  in  the  horizontal  velocity  field  as 


compared  with  the  vertical  field  and  in  part  to  the  less 
steep  vertical  concentration  gradient  above  the  peak. 

(2)  The  phase  of  the  response  relative  to  the 
gravity  wave  varies  continually  with  altitude  and  is 
180°  out  of  phase  well  below  the  peak  but  comes  in 
phase  well  above  the  peak,  in  agreement  with  the 
original  results  of  Thome  (1968)  . 

(3)  On  the  whole  (but  not  always)  the  structures 
produced  by  the  OH  emission  profile  are  more  pronounced 
below  the  peaks  of  hydrogen  and  ozone  than  above  and 
pronounced  dark  areas  can  occur  immediately  below  the 
peak  OH  intensity. 
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Figure  Captions 


Gravity  wave  model  for  a  horizontal  phase  velocity 
of  302  rn.sec."1  AW  and  AL'  are  the  vertical 
and  horizontal  velocity  fields  respectively. 


Fig.  2 


Gravity  wave  model  for  a  horizontal  phase  velocity 
of  33  m.sec"!.  A1'/  and  Ail  are  the  vertical 
and  horizontal  velocity  fields. 


Fig.  3  Unperturbed  ozone  number  density  profile. 


Fig.  4  Unperturbed  hydrogen  number  density  profile. 


Fig.  5(a)  A  comparison  of  the  variations  in  ozone  number 
density  AN  (0^)  and  total  number  density 
ANt  produced  by  a  gravity  wave  with 

a  horizontal  phase  velocity  of  302  m.sec.'l 

Fig.  5(b)  A  comparison  of  the  variations  in  hydrogen  number 
density  AN  (H)  and  total  number  density  ANj 
produced  by  a  gravity  wave  with  a  horizontal 
phase  velocity  of  302  m.sec.'1 

Fig.  6(a)  A  comparison  of  the  variations  in  ozone  AN  (0^) 
and  total  number  density  An„  produced  by  a 
gravity  wave  with  a  horizontal  phase  velocity 
of  78  m. sec . *  * 

Fig.  6(b)  A  comparison  of  the  variations  for  hydrogen  AN  (H) 
and  total  number  density  A Nj  produced  by  a 
gravity  wave  with  a  horizontal  phase  velocity 
of  78  m. sec . " * 

Fig.  7(a)  A  comparison  of  the  variations  for  ozone  Z\  N  (0,) 
and  total  number  dcnsityZX  Ny 
produced  by  a  gravity  wave  with  a  horizontal 
phase  velocity  of  60  m.sec.'l  The  two  curves 
are  out  of  phase  below  90  km.  but  begin  to  get 
in  phase  above  this  height. 


Fig.  7(b)  A  comparison  of  the  variations  for  hydrogen  AN  (H) 
and  total  number  density  A  N-j-  produced  by  a 


gravity  wave  with  a  horizontal  phase  velocity 
of  60  m.sec.'l  Again,  the  two  curves  are 
out  of  phase  below  86  km.  and  begin  to  get 
in  phase  above  this  height. 


Fig.  8  A  comparison  of  the  variations  for  ozone  An  (0  ) 
and  total  number  density  &Ny  produced  by  ^ 
a  gravity  wave  with  a  horizontal  phase 
velocity  of  33  m.sec.“l  For  this  case,  the  two 
curves  are  obviously  out  of  phase  below  90  km. 
and  are  quite  evidently  in  phase  at  high  alti¬ 
tudes  . 


Fig.  9  A  comparison  of  the  variations  in  hydrogen  A N  (H) 
and  total  number  density  A  N-r  produced  by  a 
gravity  wave  with  a  horizontal  phase  velocity 
of  33  m.sec.'^  Again,  the  two  curves  are 
clearly  out  of  phase  at  low  altitudes  and  are 
perfectly  in  phase  above  90  km. 


Fig. 10(a)  The  same  curves  as  those  in  Fig.  8  are  plotted 
using  Earl  Good's  data  for  unperturbed 
ozone  profile. 

Fig. 10(b)  The  same  curves  as  those  in  Fig.  9  are  plotted 

using  again  Earl  Good's  data  for  unperturbed 
hydrogen  profile. 

Fig.  11  OH  intensity  profile  in  the  presence  of  a  gravity 
wave  with  a  horizontal  phase  velocity  of 
33  m.sec.'l  Zimmerman's  data  are  used  for 
the  unperturbed  ozone  and  hydrogen  profiles. 


Fig.  12  OH  intensity  profile  in  the  presence  of  a  gravity 
wave  with  a  horizontal  phase  velocity  of 
33  m.sec.'l  Earl  Good's  data  are  used  for 
the  unperturbed  ozone  and  hydrogen  profiles. 
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